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Phase structure of CP^ 1 model with topological term * 
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Qpjv-i moc j e j w jth topological term is numerically studied. The topological charge distribution P(Q) is cal- 
culated and then transformed to the partition function Z(0) as a function of 9 parameter. In the strong coupling 
region, P(Q) shows a gaussian behavior, which indicates a first order phase transition at = 7T. In the weak 
coupling region, P(Q) deviates from gaussian. A bending behavior of resulting F(0) at 8 ^ n, which might be a 
signal of a first order phase transition, could be misled by large errors coming from the fourier transform of P(Q). 
Results are shown mainly for CP 3 case. 
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1. Introduction 

It is well known that CP^ 1 model shares with 
QCD many dynamical aspects such as asymptotic 
freedom, confinement, dynamical mass genera- 
tion and non-trivial topology. We are concerned 
with the topological aspects of the model and 
study the nature of the 9 vacuum. It is deeply 
associated with the non-perturbative nature of 
the strong interaction. The strong CP problem 
is one of the issues to be clarified nonperturba- 
tively. It is also expected that a new parameter 
introduced into the theory, in general, would in- 
duce a rich phase structure. Actually, in 1982, 
Cardy took, as a toy model, the Z(N) gauge model 
and by making a renormalization group argument 
showed that very rich phase structures emerged 
0. In this talk we are concerned with the dynam- 
ics of the 9 vacuum of the CP^ -1 model with 
the topological term, and report on the results 
mainly of CP 3 . CP 3 model is semiclassically free 
from dislocation. 

From the numerical point of view, the topologi- 
cal term introduces complex Boltzmann weight in 
the euclidean space time. It prevents one from ap- 
plying straightforwardly the standard algorithm 
of the Monte Carlo simulations. This problem 
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can be circumvented by Fourier-transforming the 
topological charge distribution P(Q) It is 

then necessary to calculate P(Q) in a high preci- 
sion because Fourier transformation will generate 
large error propagation. In order to achieve a high 
precision, we employ multi-histogram method as 
well as reweighting method. 



2. Formulation 



The lattice action of the CP^ 1 model is given 



by 
S = 



\-fj,Un,fj,Zn', ~r"C.C. 



(i) 



where z n is CP^ -1 variable at site n and U Uifl is 
U(l) variable sitting on a link n, fi. (3 is the cou- 
pling constant. This action is known to be supe- 
rior with respect to scaling behavior to the stan- 
dard action which is quartic in z n variables. 
We use metropolis algorithm combined with the 
overrelaxation, which is applicable to arbitrary 
N. The 9 term is added to eq.(Q) . 

S e = S - %9Q (2) 

where the topological charge is computed accord- 
ing to the definition 



vO. ( 3 ) 
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where 9 n ^ is the phase of U n ^. 

In order to avoid complex Boltzmann weights, 
we adopt the algorithm by which the partition 
function is given by the Fourier transform of the 
topological charge distribution P{Q) 



Z{9)=^e l8Q P{Q) 
Q 



(4) 



The distribution P(Q) is calculated by the real 
Boltzmann weight 



P(Q) = / [dzdz] Q e- s / / [dzdzje 



(5) 



where [dzdz]® is the constrained measure in 
which the value of the topological charge is re- 
stricted to Q, and P(Q) is normalized such that 
^2qP{Q) = 1. Error propagation in the fourier 
transform (Q) , in general, causes large errors to 
Z(9), free energy and its derivatives. In order 
to acquire reasonable results within tolerable er- 
rors, we then have to calculate P(Q) in a very 
high precision. For that purpose, a combined 
use of the multi-histogram method (set method) 
and reweighting method (trial function method) 
is made. Since P(Q) is an even function of Q, 
we divide the range of positive (and 0) values of 

Q into sets Si, i = 1, 2, 3, Each of the sets Si 

consists of 4 bins Q — 3i — 3, 3i — 2, 3i — 1, 3i so 
that adjacent sets overlap at the edge bins of each 
set, Q = 3k(k = 1,2,3,.. .). For negative values 
of Q, we use P(Q) = P(-Q). 



3. Results 

3.1. Strong coupling region 

Topological charge distribution P{Q) is calcu- 
lated for various volumes V and coupling con- 
stants. Typical behavior of P(Q) in the strong 
coupling region is plotted as a function of Q 2 in 
Fig.l for a fixed j3. Solid lines are gaussian fit- 
tings to the data, and it is clearly seen that the 
fittings work very well for all V. The coefficient 
a 2 , read off from the fit P(Q) = Aexp (— a 2 Q 2 ), 
shows that a 2 oc l/V = C/V (C is a constant), 
which turns out to be a finite size scaling law of 
the first order phase transition at 9 = tt. For that, 
the partition function is now expressed as the 
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Figure 1. P(Q) at = 0.5 for CP 3 



third elliptic function $3, Z{6) oc i?3(^,r), where 
exp(i27ri^) = e lS and exp(i7rr) = exp(— C(J3)/V). 
Infinite product 
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yield zeroes at ( = —p ^ 2n 1 \—p 2n 1 or 6 = it ± 
. (2n-i)c ( ^ — i ; 2, 3, . . . ), where ( = exp(i27ri/) 
and p = exp(«7rr). Accumulation of zeroes ac- 
cording to the finite size scaling law 1 /V indicates 
a first order phase transition at 9 — tt || . 

Resulting free energy density F(9) = In Z(9) /V 
develops a cusp at 9 = tt as V increases because 
of the periodicity F(9 + 2n) = F(9). Its behav- 
ior agrees with the result of the strong coupling 
expansion by Seiberg||l0||. 

3.2. Weak coupling region 

In the weak coupling region, P(Q) devi- 
ates from the gaussian. We tried various fits 
to lnP(Q) such as Polynomial fits, or adding 
Q 1 / 2 ,^ 3 / 2 etc. to them in order to make the 
fourier transform efficient. However, none of our 
trials was successful in the sense that resulting 
errors of F(6) is so large that meaningful results 
could not be obtained. Although there may be 
some other better fits, we decide to use the data 
itself for the Fourier transform. 

Typical behavior of the resulting free energy is 
shown in Fig. 2. F(9) develops a sharp bend at 
9b 7^ 7r. At the conference, we reported on a pos- 
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sibility of first order phase transition at 9 tt. 
Thereafter, however, we came to realize an insuf- 
ficiency of our way of analysis. So the conclusions 
we drew in the talk might have been somewhat 
misleading. The partition function receives large 
errors from P(Q), particularly from P(Q) in the 
first set Q = 0, 1, 2 and 3, because of rapidly de- 
creasing function. Z(9) becomes then meaningful 
only for the values of 9 such that Z{9) > \SZ\, 
where 5Z is a fluctuation coming from the fourier 
transform of P(Q) + SP(Q). By setting 9 = 9 fl 
such that Z(6fi) « \8Z\, data seem to indicate 
9b ~ 9fi. Since 9fi is associated with the errors 
of P(Q), its location depends on the number of 
statistics. Our preliminary analysis seems to sup- 
port this feature. 
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Figure 2. F (9) at (3=1.5 for CP 3 . Errors are not 
included. 



4. Comments 

We need to be more careful with the weak 
coupling region. Bending behavior of F(9) at 
9 =/= 7r in the weak coupling region might not lead 
straightforwardly to the conclusion that first or- 
der phase transition occurs at 9 ^ tt 0. Our 
results of volume dependence of 9b for various /3 
is quite similar to that obtained by Schierholz. 
In his papers Q Schierholz concluded that 9 C lin- 
early decreases in V, and that its extrapolation 
to V — > oo may lead 9 C — > as a possible solution 



to the strong CP problem. Our data are, how- 
ever, too noisy to draw such a conclusion, and the 
detail will be reported in |ll[] . 

In addition to CP 3 model, we reported about 
the results of CP 1 and CP 2 models with the stan- 
dard actions, quartic in z variables. In case of 
CP 1 , fixed point action is also used for compu- 
tations. Compared to the standard action, fixed 
point action shows that strong coupling behav- 
iors are seen only in the restricted region very 
close to j3 = 0. As far as the behavior of F{9) is 
concerned, these models show similar behaviors 
in both of the strong and weak coupling regions. 

The authors thank R. Burkhalter for useful 
discussions. This work is supported in part by 
the Grant-in- Aid for Scientific Research from the 
Japanese Ministry of Education, Science and Cul- 
ture (No. 10640276). 
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